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SUMMARY 


A multi-level grid method has been studied as a possible means of accelerat- 
ing convergence in relaxation calculations for transonic flows. The method 
employs a hierarchy of grids, ranging from very coarse (e.g. U x 2 mesh cells) 
to fine (e.g. 6k x 32 ); the coarser grids are used to diminish the magnitude of 
the smooth part of the residuals, hopefully with far less total work than would 
be required with, say, optimal SLOR iterations on the finest grid. To date the 
method has been applied quite successfully to the solution of the transonic 
small-disturbance equation for the velocity potential in conservation form. Non- 
lifting transonic flow past a parabolic-arc airfoil is the example studied, with 
meshes of both constant and variable step size. 


INTRODUCTION 


The multi-level grid method, for accelerating convergence in relaxation cal- 
ciilations , has been shown to be very efficient for solving elliptic problems with 
Dirichlet boimdary conditions. For background and historical material, see ref- 
erences 1 to .The idea of the method is based on the fact that in many typical 
elliptic boundary- value problems, the error is composed of a discrete spectrum of 
wave lengths , which range from the width of the region down to the width of a 
mesh cell. The short wave-length components of the error are usually diminished 
quite rapidly in a relaxation calculation, while the long wave-length components 
diminish very slowly. After only a few iterations the residual will be smooth, 
since the short wave-length error components have been eliminated; and thus the 
residual can be represented accurately on a coarser mesh. An equation called 
the "residual" equation is then solved on the coarser mesh, and the resulting 
correction is added to the last approximation on the fine mesh, yielding a sig- 
nificant improvement with very little work. 
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Since relajcation methods are currently the most attractive for obtaining 
numerical solutions to transonic aerodynamics problems , the question arises 
as to whether a multi-level, or multi-grid (MG), method can be used in a mixed 
flow with shock waves. In this paper we report some early results using the 
MG method to solve a simple transonic problem: we consider the transonic small- 

disturbance equation for the velocity potential, for nonlifting flow past a 
parabolic-arc airfoil. 


PEOBLEM DESCRIPTION 


The transonic small-disturbance equation for the velocity potential can be 
written in conservation form as : 
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Equation (l) is to be solved subject to the boundary conditions that the distur- 
bance potential, c(), vanishes at infinity and the flow is tangent to the airfoil 
surface, in the interval |x| < 1/2; i.e., 

at y = 0, (p = F’(x) for \x\ < 1/2 

^ (5) 

= 0 for, |x| > 1/2 

where F(x) is the (upper surface) thickness distribution function, T is the 
usual thickness ratio, and y? ^ ratio of specific heats, 

free-stream Mach number, and transonic similarity parameter, respectively. The 
form of equations (l) to (5) is a correctly scaled transonic similarity form, 
in that all quantities are of order 1. Equation (l) is of hyperbolic or elliptic 
type depending on whether 

U = K - (y + l)^2 (6) 

is negative or positive, respectively. 


1360 



Finite-Difference Equations 


Murman’s conservative difference scheme (ref. 5) can he conveniently pre- 
sented in terms of Jameson^s "switching function" (ref. 6) as follows: 


/l - y. A P. , + y. ^ . P. . + Q. . = 0 
[ ij 1-1, j 1-1, J ij 


(7) 


where 
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ij ij' 
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and where 


y. . = 0 if U. . > 0 
ij 10 


= 1 if U. . < 0 
10 ■ 


( 11 ) 


It should he noted here that, in' the interest of simplicity, we have pre- 
sented only the constant-step-size (unstretched grid) form of the difference 
equations. The actual computer program is written for a stretched grid, with 
the identity transformation (constant step size) included as a special case. 


Vertical Line Relaxation 

A vertical line relaxation scheme for solving equation (7) "by iteration can 
he written as : 


where 
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<f) denotes a "new” value of (p, obtained during the latest iteration sweep, 
while (J) is the value from the previous sweep. Ri j , which is the left-hand 
side of equation (7)» is evaluated with "old" values of 4>i j » as are the itera- 
tion coefficients A through E, which are given in the appendix. 


Mvilti-Grid Approach 

Residual equation .- Let us introduce a sequence of grids Gq, G2, ...» Gj^, 
where for simplicity, hj^ = 2hi[+]_, and hj^^ represents the step size of the 
Gjj grid. We can represent the iteration operator (e.g., eq. (12)) on the 
finest grid Gj,^ as: 


vhere d)., is the exact discrete solution on the grid* We can write 
M M 

♦m = “m * ''m 

where is the approximate solution and v^^ is the error. Then we have the 

residual equation ; 


= -^ 


(16) 


where is the residual of the approximation on the G^ grid. is in 

general different from in the nonlinear case, which complicates matters. 

Nevertheless, if is smooth, the error will be smooth, and the residual 

equation (l6) can be solved on a coarser grid. Thus, for example, we can write 

^M-l ^ 


£ 

where j_ is an approximation to the error v^^ on the 

denotes interpolation from the G^ to G^. After solving the problem (17) 

(usually with homogeneous boundary conditions), we interpolate the function 
1 onto the G^^ mesh, and thus form an improved approximation: 

new “ old ^M-1 


1362 


In the complete MG algorithm, the solution of equation (17) is also per- 
formed hy relaxation; and if the convergence rate falls "below a prescri"bed 
level, we can apply a similar procedure, hacking up to the Gm_ 2 grid level, 
and so on, until we arrive at G^, if necessary. The G^ grid is so coarse that 
a direct solution could he used economically , hut we have used Iteration here 
also. 


Full approximation .- In the general nonlinear case, the form of the 
operation L can he quite complicated — more so than the original operator , 

L — and thus applications to, say, the full potential equation may he tedious 
to program. It turns out that for the transonic small-disturhance equation, 
the joh is simple, and our first program did use the exact expression for L 
in an efficient way. However, there is an equivalent, easier method for solving 
the residual equation, which we call the full annroximation method, as follows: 

Suppose we add to both sides of equation (17) the f\inction 
^-1 " ^ M -1 “ ^-1 


Then, since 

^M-l ^^M-1^ ''' Sl -1 “ ^-1 

we have 

^-1 


We can now use the original operator on all the grids , which greatly simplifies 
the programming. The right-hand side of equation (20) 'is the difference 
between the residuals of u^ calculated with the coarse- and fine-grid opera- 
tors. Note that when the solution converges on the grid, then 


^ 0 (21a) 

hut will remain finite, since is a solution on the Gj^ grid; ^ 

is essentially the truncation error of the ^ operator. 
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After equation (20) is solved to sufficient accuracy, we determine the 
function 

^M-l “ 


hy subtraction at all points of the grid G , and then interpolate w to 
the grid as before in equation (l8). ~ 


RESULTS AND DISCUSSION 


In order to estimate the efficiency of the method, a work xmit can be 
defined as the amo\mt of computational effort required for one relaxation sweep 
on the ( finest ) G grid . Thus a relaxation sweep on the G, grid costs 

n = (l/^) work units, for example. Likewise, when we calculate the resid- 
uals for the Gt^ grid, we perform these calculations at the points of the G 
grid, i.e., 1/4 ^ " ^ k-x 


as few points; hence each residual calculation costs less 


than l/|+^^he effort of a relaxation sweep on the ^ grid, or approximately 
(l/4)^” • Note that this is an overestimate, since the tridiagonal system (l2) 

is not inverted, nor do we calculate the iteration coefficients during the 
residual calculations. On the other hand we did not count the work of inter- 
polation in equation (l8), for example, or any other "overhead” of that type. 

An overall estimate of efficiency can he given hy the effective spectral 
radius 




1/n 


w 


(23) 


where 


I |r^ l| = norm of after first sweep on G^ 


\\\ 


= noim of R, after n work units 
,n^II ^ w 


and 




(2l^) 
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Hence the nonn we use is the root mean square of the residual on This 

number is typically about 5 to 10 times smaller than the maximum norm in 
transonic problems. We consider an approximate solution to be converged when 


I |Ey| I < C/(no. 


of grid points) 


( 25 ) 


where the prescribed constant C is typically chosen as 1 so as to estimate the 
nominal truncation error. 


Unstretched Grids 

In the case of a grid with constant steps in both directions , the present 
MG method performed quite well. Some typical results are summarized in table I 
and discussed briefly in the following. 

In all cases, the MG runs were made with a relaxation factor OJ = 1.0 on 
all grids. 

Laplace *s equation with smooth boundary conditions .- To illustrate just 
how fast the MG method works for a nice, smooth, elliptic problem, we present 
in table I results for the solution of Laplace’s equation with the prescribed 
normal derivative equal to sin ttx along y = 0. Because of the smoothness of 
the boundary data, it could be expected that interpolating a converged G^ 

(32 X i 6 grid) solution onto G^ will give a very good starting approximation 
for G 5 . This is true, for although the convergence rate on G^ yielded 
a = . 583 , the efficiency of the two combined levels is more like a = .46! In 
contrast, successive line overrelaxation (SLOR) achieved a = .924 on G^ , 
starting from the zero solution, and using a relaxation factor = 1.85 • 

Nonlinear airfoil flows .- The next three entries in table I show the 
results for the nonlinear problem of flow over a parabolic-arc airfoil. In 
these cases, the Neumann boundaiy condition is an ”N-wave” — far from smooth.— 
but the M^ = 0.7 subcritical case (i.e., no supersonic flow) converged as 
well as the previous smooth problem; hence, it can be concluded that discontinu- 
ous boundary conditions do not deteriorate MG performance. The "combined" mode 
of operation, where the converged solution for Gi| is used to start G^ , was 
not helpful, since the truncation errors around boundary singularities and 
shock waves were so large. That is, the Gi;-solution gives a large residual* 
when interpolated on the G^-mesh. The = O .85 (moderately supercritical) 
case had 124 supersonic points out of a total of 2l45 mesh points on G 5 , or 6%. 
The relative efficiency between MG and SLOR is still unaffected. In both of 
the aforementioned nonlinear cases, the SLOR runs were carried out with co = 1.85s 
which was found to be near optimal by experiment. 

The last of the unstretched grid cases is = 0.95 (highly supercritical), 
with 355 supersonic points. The flow pattern exhibited a weak oblique shock at 
the trailing edge, followed by a triangular region of nearly constant supersonic 
flow, which was terminated by a normal shock in the wake. The final number of 
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supersonic points was established after 38 work units , and the solution con- 
verged after 67.6 work units, giving a = O. 858 . The SLOR run was unstable 
with 03 = 1 . 85 , and had to be "babied" by slowly increasing 03 , using an inter- 
active remote terminal. The best result achieved was n^ = 228, with a = .957* 


Stretched Grid? 

An attractive way to satisfy the boundary condition at infinity is to 
transform the independent variables such that the infinite space is mapped 
onto a finite domain. However, it became quickly evident that vertical line 
relaxation alone is not the best way to relax the solution for a stretched grid, 
either in the MG mode or simple SLOR. Analysis of the difficulty shows that all 
the high-frequency error modes are not rapidly damped if the mesh aspect ratio 
differs significantly from 1.0; the success of the MG method, of course, hinges 
on this feature. The analysis, not given here, also indicates that a solution 
to this problem is to sweep in all directions alternately (forward, backward, 
up, and down). 

The last entry in table I shows the resxilts of a stretched-grid case, 
again for M^ = 0.95* The deterioration of the MG method is clear; some bene- 
fit over SLOR is achieved, however, by the MG method. 


CONCLUDING REMARKS 


The MG method for accelerating relaxation calcxxlations has proved to be 
applicable to nonlifting transonic flows with embedded shock waves. The method 
appears to work from three to five times faster than optimal SLOR on unstretched 
grids of moderate size (6U x 32); the relative advantage of MG over SLOR 
increases as the grid gets finer, since the MG convergence rate is nearly inde- 
pendent of mesh size. It is probable that the gains in three-dimensional cal- 
culations would be even more impressive, since each coarser grid requires 
only 1/8 the work of the next finer grid. 

On stretched grids , the present MG approach slows down , being only about 
twice as fast as SLOR. It is felt that a remedy is the use of alternating- 
direction relaxation sweeps. 

In the futiore we hope to develop the MG method for flows with lift; for 
otherwise it will have limited usefulness in aerodynamics. 

During the course of ovir work. Professor Antony Jameson of the Courant 
Institute of Mathematical Sciences, New York University, also carried out 
research on the multi-grid method. He showed independently that the "full 
approximation" approach would work, and some of his attempts at alternating- 
direction sweeps have been encouraging. Our many discussions have been 
beneficial. 
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APPEiroiX 


ITERATION COEFFICIENTS 

We have used various choices for iteration coefficients in equation (12). 
The coefficients used to make the calcvilations presented in this paper are 
simply based on the Newton linearization of equations (7)» (8)> and (lO). They 
are as follows : 

First define; (dropping the j index, since all quantities are 

evaluated at the same j ) 


'>i.i = I' - - ♦i) 

2 Ax 


(Al) 


Then we have 


2 2 


U.Ax"^ 

1 


(A2) 


A = C = -Ay"^ (A3) 

B = 2Ay“^ + 2 (1-y^) b^_l (aU) 

2 


D = (l-y^) - 2M._1 U^_l (A5) 

2 

E = y. T b. 3 (A6) 

*^ 1-1 1 -^ 

2 

where 

y. = 0 if U. > 0 
1 1 

(A7) 

= 1 if < 0 
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TABLE I.- SUMMARY OF MULTI-GRID RESULTS, 6 k x 32 CELLS 


Problem description 

Si 

Effective spectral radius for — 

MG 

SLOR 

Unstretched 

Laplace's equation^ smooth 

0.583 

0.924 

grid 

boundary conditions 

(0.U6 combined levels) 



Parabolic airfoil, = O.7O 

.5U9 

.868 


Parabolic airfoil, 

.593 

.855 


Parabolic airfoil, = .95 

.858 

• 957 

Stretched 

Parabolic airfoil, = 0.95 

0.936 

0.97k 

grid 





^See equation ( 23 ). 
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